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We propose a model of a spinless electrical charge as a self-consistent field configuration of the 
electromagnetic (EM) field interacting with a physical vacuum effectively described by the logarith- 
mic quantum Bose liquid. We show that, in contrast to the EM field propagating in a trivial vacuum, 
a regular solution does exist, and both its mass and spatial extent emerge naturally from dynamics. 
It is demonstrated that the charge and energy density distribution acquire Gaussian-like form. The 
solution in the logarithmic model is stable and energetically favourable, unlike that obtained in a 
model with a quartic (Higgs-like) potential. 
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1. INTRODUCTION 

Two of the oldest actual problems in fundamental par- 
ticle physics relate to finite self-energy and the possible 
extendedness of electrically charged elementary particles. 
This research direction is complicated by the fact that no 
experimental evidence of either the internal structure or 
spatial extent of, e.g., the electron has been found down 
to 10~^^ cm. The mere postulate that a certain amount 
of matter with mass, charge and spin can be located in- 
side a set of zero spatial measure looks implausible to a 
physicist's mind. This assumption, however, might be 
one of the reasons why unphysical divergences appear 
in quantum field theory (QFT). This difRculty already 
arises at the classical level: according to the standard 
theory of electromagnetism, the electrical field of a point 
charge in completely empty space (trivial vacuum) is de- 
scribed by the inverse square law, therefore the energy 
density of the electrical field integrated over the whole 
space turns out to be infinite. As a result, the total 
mass-energy of the point charge, together with its field, 
becomes infinite, therefore, such a system would be im- 
possible to move. At the quantum level, this problem 
manifests itself in ultraviolet divergences appearing in 
loop diagrams. In some theories, these divergences can 
be removed by means of regularization and renormaliza- 
tion procedures. This can be very useful for doing spe- 
cific computations, but does not shed much light upon 
the essence of the problem [iHH. Theoretical attempts 
towards better understanding should not, however, be 
abandoned. 

Historically the first effort to address this problem was 
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probably a model which described the electron as a ball 
with spatially distributed electrical charge. That model 
conflicted with relativity, however, because the ball was 
assumed to be absolutely rigid. The description of spin 
was also not clear. Similar difficulties were found in mod- 
els proposed by Abraham and Lorentz although 
work in that direction continues [l^. Dirac's shell 
model of the extended electron, together with its subse- 
quent modifications and variations [ll| - [l^ . provides an- 
other notable research direction. Yet another model of 
an extended electrical charge was the Einstein's worm- 
hole approach. In his model of an electron the electrical 
flux lines enter one side of the wormhole and exit from 
another, resulting in the front side looking like a neg- 
ative charge and the rear like a positive charge. The 
wormhole models were criticized by Wheeler for issues of 
stability, non-quantized charge, wrong mass-charge ra- 
tio and spin [1^. Numerous attempts towards flnding 
regular particlc-likc solutions were made in conventional 
and nonlinear electrodynamics [l6l. both with and with- 
out engaging general relativity jl7H28| . Another inter- 
esting approach is wave-corpuscle mechanics (for a re- 
view see [29| ) . The general mathematical formalism used 
therein formally resembles the one used in this paper. 
The important difference, however, is about underlying 
physics: the origin and explicit form of their nonlinear 
self-interacting wave term Gd'f/'p) are not specified on 
physical grounds and a fully satisfactory particle-like so- 
lution is not given. 

A popular approach to the classical electron model was 
made using the Einstein-Dirac [isj and Kcrr-Newman 
(KN) solutions [30l - l42j (an extensive bibliography can be 
found in [4^ ) . While the original KN solution does have 
the correct gyromagnetic ratio, it also contains a naked 
singularity and thus requires an additional mechanism 
to circumvent the regularization problem; the story is 
far from being complete yet. 
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Intuitively, from a fundamental theory one would ex- 
pect that spatial extent is not ab initio built-in, but nat- 
urally emerges from dynamics. In this paper we propose 
a model of a charged particle whose spatial extent, ob- 
servable charge and mass emerge as a result of the in- 
teraction of the EM field with the physical vacuum. For 
simplicity we neglect internal degrees of freedom, such as 
spin, isospin, etc., so that wc can assume spherical sym- 
metry where possible. The resulting solution describes 
a charged object which does not have a boundary in a 
classical sense; its stability is supported not by surface 
tension but by nonlinear quantum effects in the bulk. 
This makes our model more realistic from the quantum- 
mechanical point of view, since the actual observability of 
a definite boundary with smooth surface would be as con- 
tradictory to the quantum uncertainty principle as the 
notion of a smooth trajectory or worldline in the quan- 
tum realm. This can be shown by performing a simple 
Gedankenexperiment: making an extended object with a 
definite boundary propagate through space and measur- 
ing the velocity and position uncertainties on its surface. 
In turn, it means that the surface tension is a well-defined 
notion only in the classical limit, but for more fundamen- 
tal purposes it must be used with utmost care. 

The structure of the paper is as follows. The phc- 
nomenological approach to physical vacuum is described 
in the next section, the main equations of the model can 
be found in Section [31 the regular solution and its prop- 
erties are analyzed, both analytically and numerically, 
in Section |4l A comparison between our solution and 
its Higgs-type (quartic) counterpart is done in Section [Sj 
and conclusions are drawn in Section [51 



2. PHYSICAL VACUUM 

As mentioned above, the Coulomb divergence problem 
essentially means that one cannot find regular particle- 
like solutions of the Maxwell field in empty space, not 
even in general relativity From a quantum physi- 

cist's point of view, however, this problem is not as se- 
vere as it looks to a non-quantum theorist, because the 
notion of absolutely empty space (or "mathematical vac- 
uum") cannot be realized in nature anyway. This is be- 
cause the existence of such space seriously contradicts 
quantum-mechanical laws. According to the latter, the 
genuine (physical) vacuum must be a non-trivial quan- 
tum medium which acts as a non-removable background 
and affects particles propagating through P, . 

At this time, no commonly accepted theory of physical 
vacuum exists. The amount of experimental and obser- 
vational data is still too far from conclusively identifying 
a single model. One of the candidate theories lies within 
the framework of the superfluid vacuum approach [46l - 
[Sll . This theory is based on the idea [52| - [53 | that a phys- 
ical vacuum can be viewed as some sort of background 
superfluid condensate described by the logarithmic wave 
equation (the latter was studied previously on grounds 



of the dilatation covariance [s^ or separability (s^-fs^). 
It was shown that small fluctuations of the logarithmic 
condensate obey the Lorentz symmetry and can be inter- 
preted not only as the relativistic particle-like states but 
also as the gravitational ones [sJl depending on a type of 
mode. In this approach, therefore, the Lorentz symmetry 
is not an exact symmetry of nature, but rather pertains 
to small fluctuations of the physical vacuum, and thus 
gets deformed at high energies and/or momenta. 

As long as the superfluid vacuum approach must be 
fully consistent and applicable to reality, one would ex- 
pect that if the empty space is replaced by the logarith- 
mic vacuum condensate then the behaviour of the con- 
ventional Maxwell field becomes regular in the presence 
of particle-like solutions. This behaviour is going to be 
the main subject of the current study. One would also 
expect that the spatial extent mentioned above must ap- 
pear naturally in the approach. This has already been 
shown at the non-relativistic level in (soj . so in the cur- 
rent study we will investigate the relativistic case. 

The effective low-energy Lagrangian for the Maxwell 
field interacting with small fiuctuations of the physical 
vacuum was proposed in [s^ , along with a mass genera- 
tion mechanism which was analogous to the Higgs one. In 
this paper we propose a different mass generation mech- 
anism which uses the same Lagrangian, except that the 
scalar potential is assumed to be "upended" . As will be 
shown below, the latter solves the issue of a wrong sign 
in the quadratic (mass) term of the potential at energies 
above the symmetry-breaking scale - when symmetry is 
unbroken and false vacuum is stable. In the case of three 
spatial dimensions the action is proportional to / d'^x C 
where, adopting the natural units c = fi. = eo = 1 and 
metric signature (H ), we assume 

C = -^F,^F^^+^^\D,^f~Vm'), (1) 

where = 9^ + ig^^, Ff^^ = 9^^^, - a^A^, and the 
vacuum-induced field potential is defined as (up to an 
additive constant) 

Vm') = [ln(a3|*n - 1] + a-'} , (2) 

where d and /3 are parameters of dimensionality length 
in adopted units. In the underlying theory of superfluid 
vacuum the parameter /3 is related to the quantum (non- 
thermal) temperature which is conjugated to the Everett- 
Hirschmann information entropy |59j whereas a can be 
related to the characteristic inhomogeneity scale of the 
superfluid j60i|. The potential ^ is regular in the ori- 
gin - while the logarithm itself diverges there, the factor 
l^'P recovers regularity. There is a local maximum at 
1^*1 max = d~^^^, i.e., it always has the (upside-down) 
Mexican-hat shape if plotted as a function of see Fig. 
[T] In what follows we call this potential logarithmic - due 
to the property dV/d\'^\'^ oc ln(a^|*|2) which yields the 
logarithmic term in the corresponding field equation. 

Wc emphasize that, according to this approach, the 
Lagrangian ([T]) is an approximate one, thus it is not valid 
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FIG. 1; The field-theoretical potential ((2| (in units of /3a^) 
versus Re(v[') (in units of a^''''^). Vertical lines condition- 
ally represent inequality (|3}- In an approximation when the 
symmetry-breaking energy scale is much less than the vacuum 
one (> 10 TeV) these walls can be assumed infinite. 



for arbitrarily large (or short) scales of energy (or length) 
which means that the 4* field cannot take arbitrarily large 
values: 



1*1 ^ l^-J < oo, 



(3) 



where j^'d = lim |4'| is some limit value corresponding 

E^Ea 

to the cutoff energy scale Eq which is also the character- 
istic energy scale of the vacuum. In the effective theory 
the appearance of upper bound for will be shown in 
Section 14.21 below, whereas in a full theory it could re- 
sult from, for instance, the normalization condition for 
^, similarly to the one in the theory of Bose-Einstein 
condensation. The constraint ([3]) also means that the 
potential ([2]) does not have to be bounded from below, 
as in a standard relativistic QFT. Alternatively, one can 
take the potential ^ with an opposite sign, thus mak- 
ing it bounded from below at positive /?, but treat * as 
a phantom field. 

Further, performing the rescaling 



= V^*, ^ = a/3, a = 0^/3, (4) 
we can rewrite ([T]) and ([2]) in a more regular form; 



(5) 



m = -r' {l^P [inla^lV^n - 1] + ^} . (6) 

Note that by expanding the potential in the vicinity of 
= e/a^, £ being a non- negative dimensionless num- 
ber, one arrives at the following perturbative expression 
(up to an additive constant): 



F(^) 



Aefi 



1 



where Acff = —12a^/el3 is the effective quartic coupling 
and nit = \/2{l — lne)//3. If the radicand of nib is non- 
negative then TOb can be interpreted as the mass of an 
effective scalar particle (before the symmetry breaking). 
Indeed, one can always quantize the approximate model 
by analogy with a quartic (hence renormalizable) scalar 
QFT in the vicinity of a non-trivial vacuum represented 
by the ground-state solution of the original model (|6]) 
which will be discussed in the following sections. To date, 
a number of different quantization approaches have been 
developed for such cases - see, e.g., works j6ll463j and 
references therein. 

We have now expressed the physical parameters of the 
scalar sector, such as mass and coupling, in terms of the 
primary parameters of our theory. If the value of e is 
close to one (or, at least, less than the base of natural 
logarithm) then it is indeed important that the potential 
dS]) has the upside-down Mexican hat shape (cf. a > 
and /3 > 0) otherwise the quadratic term would appear 
with a wrong sign. Also the effective quartic coupling 
turns out to be negative in this case which is a remark- 
able difference from the standard Higgs potential and 
thus it can serve for experimental testing. In any case, 
the interpretation based on ([7]) is only approximate (for 
instance, such series expansion does not converge to (|6]) 
for very small {ipD, therefore, in what follows we will be 
working with the exact expression for V. 



3. FIELD EQUATIONS 

The field equations corresponding to the Lagrangian 
([5]) are given by 



(8) 



dV 



D^D^^i^ +Q^= [D^D^^ + In (a^lV'H] ^ = 0, (9) 
where the current j'^ is defined as 

f = ^g[{D'^^Pr^P-r{D^i^)]■ (10) 
We look for a solution in the electrostatic form 



A,, 



0(r),O , V^(t,r) = c-'^*^(r), (11) 



with E being a real- valued constant. Then the equations 
of motion become simply 



ip" + -^p' = -(£:-g0)^V-^"Vln(a^V'^), (12) 
r 



+ -0' = ~2g{E~g^)^\ 
r 



(13) 
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where primes indicate derivatives with respect to r. In- 
troducing the quantities 



-g<P\/^ = -#v a/3, 



E = = eJ&P, ~g = gJ ^ - 



9\ 

a 



we obtain 



X 

2 ~, 



_\ 2 _ 



2.92 {e+~a 



(14) 
(15) 

(16) 
(17) 



where primes indicate derivatives with respect to x when 
apphed to tilded quantities. One can see that equations 
depend on only one parameter g, whereas E = E{g) can 
be treated as an eigenvalue at a given g. The full the- 
ory of physical vacuum would provide the value of E (or, 
equivalently, "0(0)1 see the analytical solution section be- 
low) as a function of the primary parameters (from, e.g.. 



some sort of normalization condition, cf. (S^) but in the 
approximate theory E stays a free parameter which can 
be fixed only from external considerations. 

Further, for a given solution, the energy density e is 
defined as 

e = iDo^*Z?V + \D,rD'i^ + V{^) + \\£\\ (18) 

where i = 1,2,3 and £ = — is the electric field 
strength. Substituting the ansatz (fTTj) wc obtain 



l3a'e^-[^E + (b) ^' + -0' - V' In (^^^ j + ^ v • 

^ (19) 

00 

The total energy can be calculated asW = ATT J r'^e{r)dr 



so we obtain 



where we denoted the dimensionless total energy 



(20) 



E + 4>Y tp^ + i^'^ 



202 In (V-^) 



X dx. 



(21) 



r 



For an observer in the reference frame associated with 
the center of mass of a localized solution the quantity 
W/c^ is equivalent to the rest mast of a corresponding 
particle. 



4. PARTICLE-LIKE SOLUTION AND ITS 
PROPERTIES 



4.1. Asymptotic behaviour 

To search for the regular solution, the functions V'(x) 
and (j^ix) should have the following behaviour near the 
origin: 



When substituting this into (fTB]) and (fTT]) we obtain the 
solution 



7 00 
^2 = -y 



62 = -g^ [E + 4>o) i'l 



(24) 
(25) 



Further, the asymptotic behaviour at a; —> 00 is given by 
^ ^ e^(3-s=-=) [1 + o{g')] , ^ (26) 



where q is some constant to be determined. It is instruc- 
tive to relate the bare charge g to the observable one 
q = q/g. By integrating ([T^ we obtain 



r 

'£ = 2gJ iE-g^)i;\^dr, 



(27) 



^(x) = ^p„+i,,- + 0{x''), (22) 

2 

0(x) = 0o + 02^ + O(x4). (23) 



and in the limit r — > 00 we have r^£ — > q hence 



q = 2g 



00 

/ (E-gcj,) 



ijp'r^dr. 



(28) 
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Thus, we arrive at the foUowing relation between the bare 
and observable charges 

q = ^iCg) = '^-g^l^m, (29) 

where we denoted 

oo 

I{g) = j {e + ^^ iP-dx. (30) 



Expression ([26|) shows us that at large distance we re- 
cover the Coulomb potential while the field V' decreases 
exponentially. Thus, the field ■0 is in fact unobservable, 
unless very short length scales are probed. From the 
asymptotics of the solution one can infer that the charge 
radius of the solution is determined by the parameter /?: 

size - - (31) 

which essentially means that the combination of param- 
eters /? = a/3 must have an extremely small value for the 
known elementary particles. For instance, if one takes 
the values of the classical radius e^/m as conservative es- 
timates, then for the electron and muon one would obtain 
constraints of /3(e) < 10~^^ cm^ and /3(^) < 10^'^^ cm^, 
although it is not entirely clear whether the classical ra- 
dius should be analogous to the "smearing" size which is 
our definition of size here. 



4.2. Approximate analytical solution 

While the exact expression for a full analytical solution 
is unknown, it is possible to solve the system ([T6|) and 
(|17p using the approximation of weak EM coupling 

f « 1, (32) 

which is equivalent to ^ a^//3 or <^ a//3. The ob- 
servational constraints suggest that this approximation 
might have a good chance to be valid for the known el- 
ementary particles - unless a turns out to be very small. 
On the other hand, as long as our approach is an effec- 
tive one it has certain applicability conditions - and one 
of them is that the vacuum effects predominate the elec- 
tromagnetic ones. Therefore, large values of g might push 
our approach outside its applicability range and thus the 
corresponding approximation is not very interesting from 
the physical point of view. 

Thus, imposing the boundary conditions 

<^(0) < CX3, i{+oo) 0, ^^(-l-oo) < oo, i>'{{)) = 0, (33) 

one obtains a solution which is regular for ^ r ^ +oo 
(see appendix for the details of derivation) : 



1 -f ^ (2x2 + 1) e""erf(a;) 
2x 



org'), 



where dimensionless energy E can be also expressed via the boundary value 

1 



(34) 
(35) 

(36) 



with the square root being defined up to a sign. Of course, this formula is valid only if the magnitude of is bounded 
from above: 



1^1 l^ol 6^/2, 



(37) 



which a posteriori affirms the condition of applicability (|3]), although this upper bound does not necessarily saturate 
the critical value there: j^d ^ (e/a)^/^. 

Further, one can check that the electric part indeed has the Coulomb behaviour at large r, and then the effective 
charge can be computed as 



therefore, the observable charge. 



q = q/g'' 



q^\^fE^^-^' +0{g^ 



(38) 



(39) 
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depends on the whole combination of parameters describing the interaction of the electromagnetic field with a physical 
vacuum. The dimensionless total energy (j2ip turns out to be 



W 



_3_ 

16' 



'IT e 



2E^ 



1 



3^2 



E- 



(l2E'^ - is) 



0{g 



(40) 



which can also be written in terms of the observable 
charge q and rest mass W: 



W kW, 



(0) 



l + 4^/2/3-laVe^'-^^i:^l, (41) 



where 



Wi 



(0) 



S TT^/^V? 



2 a3 V^"2'^^ 



(42) 



so one can see that the obtained formula does not contain 
any divergences. 

It is also apparent that mass W does not vanish when 
charge is set to zero, which indicates that the theory is 
also capable of incorporating non-charged particles into 
the scheme, by taking the corresponding limit. In fact, 
the mass formula (pij) implies that for an electrically 
charged particle with mass W there can exist not only 
an antiparticle of the same mass but also a neutral par- 
ticle of related mass W(o) • It is interesting that the ratio 
VF/VF(o) grows exponentially with growing \E\, which re- 
sults in two possible scenarios: (i) the mass of a neutral 
partner is very small (yet non-zero) as compared to the 
mass of a charged one: this happens if \E\ ^ 1 (or, 
equivalently, \tpo\ <C 1); (ii) if \E\ is of order one or less 
then both masses would be of the same order of magni- 
tude. The possible phenomenological implications of this 
mechanism are discussed in the conclusion. 



0.00 - 




FIG. 2: Profiles of ^{x) (top curve) and tlie electrostatic po- 
tential 0(2;) (bottom curve), computed for eigenvalue E — 
2.8436935588. 



4.3. Numerical solution and stability 

While we have managed to find the approximate reg- 
ular solution analytically, it is important to check that a 
regular solution exists for non-small ^'s and that terms 
with higher-order powers of g will not introduce any spa- 
tial singularities. For this purpose we solve equations 
pop and (|17p numerically. For the computations we 
choose 5 = 1 and the following boundary conditions 

1^(0) = 0.1, 7/''(0)=0, 0(0) = -0.1, 0'(O) = O, 

(43) 

and E is treated as an eigenvalue. It should be noted 
that (t>{x) must be always taken as non-positive on the 
positive semi-axis of x, due to the asymptotic require- 
ments (pS)) . The numerical solution is presented in Fig. 
[21 and in Fig. [3] the corresponding profiles of the electric 
field £ = ~d(t>/dx and x^£ are given. From these one sees 
that the electric field is regular at the origin and asymp- 
totically displays Coulomb behaviour. The profile of the 
dimensionless energy density is shown in Fig. 21 




FIG. 3: Profiles of the electric field £{x) (top) and x^£{x) 
(bottom curve). 

The direct stability analysis of the solution is compli- 
cated by the fact that the perturbed electric field be- 
comes time-dependent, which leads to the appearance 
of a magnetic field such that this system cannot be re- 
garded as spherically symmetric anymore. It is, however, 
still possible to use the energy-based arguments as well as 
to investigate the behaviour of an effective Schroedinger 
equation potential. Let us consider the dimensionless to- 
tal energy W given by (PT|) , as well as the energy of the 
field alone, W^, = W^l^^o, and the energy of the elec- 
tric field alone, Wcf, = VF|^_>o. Then we can define the 
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FIG. 4: Profile of dimensionless energy density I3a^e(x), using 
the same E eigenvalue as in previous figures. 



FIG. 5: Efi^ective potential (|47p versus x at <?=!, using the 
same E eigenvalue as in previous figures. 



dimensionless binding energy as 



AW = W - (W^ + W^) 



U + ip'^x'^dx 



(44) 

where all the potentials are assumed to be given by our 
regular solution. One considers the following two possi- 
bilities. If binding energy is positive then it is necessary 
to add a certain amount of energy to create a regular elec- 
tric charge when coupling to the tp field. In the opposite 
case binding energy gets released during the process, and 
the energy of the whole configuration is smaller than the 
sum of the energies of the separate electric and scalar 
fields. 

Evaluating the binding energy on approximate solution 
([34|) . ([35]) . we find that it is negative-definite 



AW 



gEc 



(45) 



which means that the creation of the regular elec- 
tric charge in the logarithmic model is energetically 
favourable. 

Another way to study the stability of the solution is 
to write it as a solution of the Schrodinger equation for 
a fictitious particle 



(46) 



where A = d? /dx^^ e ^ E"^ and the effective potential is 
derived as 

Feff(x) = -2E^{x) - 4'\x) - In [i>^x)], (47) 

where the tilded potentials are given by our regular so- 
lution. According to (pS)) . the asymptotic behaviour of 
the solution implies that the effective potential Ves oc a;^ 
at large x, (see also Fig. O, and thus the "particle" is 
always localized in a finite region of x. With respect to 
the solution itself this means that it cannot spread or be 
destroyed when subjected to small perturbations. 
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FIG. 6: The profiles of the dimensionless energy W{g) (T- 
curve) and the eigenvalues E{g) (■ - curve). 



5. LOGARITHMIC VERSUS QUARTIC 
POTENTIAL 



One may wonder whether the singularity-free solution 
exists when the scalar sector of our model is controlled 
not by the logarithmic potential ^ but by a more or- 
thodox one, such as the Higgs-type (quartic) potential: 



(48) 



Corresponding dimensionless equations with the 
ansatz (fTTj) are 



X 

2 r, 



E 



- xr + 1 



(49) 
(50) 
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1.0 



at large x is given by 



FIG. 7: The profiles of ^(a;) (top) and 4'{x) (bottom) for 
the quartic model. Eigenvalue ^{Q) = 1.12345, g — 0.1 and 
parameters E = 0.1, A = 1.0 have been used. 



where we introduced following dimensionless quantities 



X = mr, ip 



iIj ~ _ E 
VKO)' m' 



m 



(51) 



The boundary conditions are 

?A'(0) = 0; 0(0) = -1.04; (^'(0) = 0. (52) 

As in the logarithmic model, the regular solution exists 
only if the potential (|48)) opens down, i.e., when >f > 0. 
The profiles of ij^ix) and (f){x) are presented in Fig. [T] For 
technical reasons in this case an eigenvalue is V'(O) not E. 
In Fig. [8] the profile of the electric field S is shown. In 
order to show that the electric field asymptotically has 
Coulomb behaviour we also present the profile x'^E(x) in 
Fig. IS] From these figures one can see that the qualitative 
behaviour of the potential 4){x) and the electric field £{x) 
are the same as for the logarithmic potential. 




FIG. 8: Electric field versus x for the quartic model. The top 
curve is lQ£{x), the bottom one is x^£{x). 

The asymptotic behaviour for the functions i}} and (p 



'tp{x) — > ip, 
4>{x) - 



(53) 
(54) 



where "000 is some constant. This still looks very similar 
to what we had earlier in the logarithmic case, however if 
we study the stability of this solution then differences do 
arise. At first, if one computes the binding energy sim- 
ilarly to pij) then it turns out to be positive, therefore, 
the creation of the regular electric charge by coupling 
electrical field to the quartic scalar one is energetically 
unfavourable. Further, if one computes the fictitious- 
particle potential for this solution, cf. (|47)) . 



(55) 



V,ff{x) = -2E4>ix) - ^^(x) - X^^ix) + 1, 



then one finds that it approaches a constant at large x 
(see also Fig. The "particle" is not, therefore, nec- 
essarily localized in a finite region of x. With respect to 
the solution itself, this means that the latter can spread 
or become unstable against small perturbations. 




FIG. 9: Effective potential (|55p versus x at g 
quartic model. 



1 for the 



6. CONCLUSION 

The classical model of a spinless electrical charge is 
described as a self-consistent configuration of the EM 
field interacting with fluctuations of a nontrivial phys- 
ical vacuum effectively represented by the logarithmic 
Bose-Einstein condensate. We have shown that a reg- 
ular solution does exist - as opposed to the case of the 
EM field propagating in absolutely empty space. In 
this regard we recall the state of affairs in quantum me- 
chanics: the Dirac/Schrodinger equation without any ex- 
ternal potential has the de Broglie wave solution; the 
Dirac/Schrodinger equation with the external electro- 
static field yields the regular wave functions (the hydro- 
gen atom being an example), but the Dirac/Schrodinger 
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equation coupled to Maxwell equations does not lead to 
a regular stationary solution. The reason is that the 
Dirac/Schrodinger equation ab initio describes a point- 
like particle which might be a good approximation for 
long- wavelength measurements, but in higher-energy and 
shorter-length regimes this approximation eventually be- 
comes too crude, since it neglects internal structure and 
non-zero spatial extent. Among other things, this leads 
to the densities of energy and charge becoming infinite at 
the particle's position. Here we have shown that by in- 
troducing an additional player on scene, the physical vac- 
uum condensate, one can obtain a regular solution, thus 
endowing particles with internal structure and spatial ex- 
tent. The solution turns out to be stable and energeti- 
cally favourable. Using its features, some observational 
constraints for the parameters of the theory have been 
derived. We also specified the conditions under which 
our model can be (approximately) interpreted in terms 
of a scalar particle and those under which it cannot. 

Further, we have established, both numerically and an- 
alytically, that the mass and spatial extent of a charged 
particle emerge due to interaction of the EM field with 
vacuum. It has been demonstrated that the average 
charge radius becomes non-zero, and the charge density 
acquires a Gaussian-like form. Looking at the form of the 
analytical solution from Section [4.2[ one can infer that it 
describes the object without border in a classical sense, 
therefore, its stability is supported not by surface tension 
but by nonlinear quantum effects in the bulk, similarly 
to the non-relativistic case [s^. Due to non-singular be- 
haviour of the solution at the origin, the derivation of 
self-energy turned out to be entirely divergence-free. 

The derived mass formula (|1T|) suggested that for an 
electrically charged massive elementary particle there ex- 
ists not only an antiparticle of the same mass (in the 
leading-order approximation with respect to the Planck 
constant, at least) but also a neutral particle of related 
mass. This might explain, at least qualitatively, why an 
electrically charged elementary particle is very often ac- 
companied by a single neutral particle of a similar kind, 
but not vice versa. Indeed, such "mass pairing" feature 
has been observed (provided one disregards the influence 
of internal degrees of freedom such as spin, isospin, etc.) 
not only for elementary particles such as leptons and 
weak bosons, but also for stable composite ones such as 
nucleons (quarks might not fit this scheme since they are 
confined inside hadrons). We presented some arguments 
for why the rest mass of a neutral partner can sometimes 
be so much smaller, yet still non- vanishing, than the mass 
of the charged one (leptons), and sometimes they are of 
the same order of magnitude (weak bosons or nucleons). 

Finally, we have compared the logarithmic vacuum 
model with one based on a Higgs-type (inverted quar- 
tic) potential. It turns out that the corresponding reg- 
ular solution is unstable and energetically unfavourable, 
in contrast with the logarithmic case. 
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Appendix: Derivation of approximate solution 

Here we provide more details regarding derivation of 
the approximate analytical solution from Section r4.2l As- 
suming ([5^ and ([55)1 . we will look for a solution of the 
system (|16p and (|17l) in series form 

^ = ^'o(x)(l+e*i(x)+0(e')), (A.l) 

where £, = g^, and $i(x) and ^^(a;) are functions to be 
determined. Keeping terms of the order C'(^) and below, 
equations (|16p and pT)) can be reduced to the following 
set of four differential equations: 

_( 

X 

^ + $0 ; ' 



1.^' + = 0, 



2^, 



+ 2 ( i + S ) ^[ + 

-f-2(£; $o)$i +2*1 =0. 
By solving them we obtain 

X 



= 0, 



(A.2) 
*o = 0, (A.3) 
(A.4) 

(A.5) 



A + e C2 



^^.[E + A)e3-(^+A)%rf(x) + O(C'), (A.6) 



2x 



^i = i(^ + A)V-(^+^)'-^' 
4 

l-h^(2a;2-t-l)c^'erf(a;) 
2x 



(A.7) 



-C2(i^ + A)+C3 [x 



1 

2x 



+ C4 



-^(2x2-l)erfi(x)-e" 



(A., 
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where Ci and A are integration constants whose values 
must be fixed by means of the boundary conditions. 
Imposing (|33|) . one obtains; ci — = 04 ^ and 
C2 = — A/^. Then, after making the energy redefinition 



E + A ^ E (alternatively, one can set A = C2 = from 
the beginning) , one eventually arrives at expressions p4p 
and ([35|l . 
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